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Abstract 

After having closely re-examined the notion of a Levy's stable vector, 
it is shown that the notion of a stable multivariate distribution is more 
general than previously defined. Indeed, a more intrinsic vector definition 
is obtained with the help of non isotropic dilations and a related notion 
of generalized scale. In this framework, the components of a stable vector 
may not only have distinct Levy's stability indices a's, but the latter may 
depend on its norm. Indeed, we demonstrate that the Levy's stability 
index of a vector rather correspond to a linear application than to a scalar, 
and we show that the former should satisfy a simple spectral property. 



1 Introduction and notation 

In order to define Levy stable (or a-stable) vectors valued on a real Hilbert 
space H and their corresponding multivariate distribution in a more gen- 
eral way than the classical and nearly isotropic definition Jl) , we will use a 
notion of scale |] x || which, as discussed in the following section, is weaker 
than the canonical norm | x | defined by the scalar product (x, y); 

l£hfe*) 1/2 (1) 
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Let L(H, H) denote the set of (continuous) endomorphisms of H. Al- 
though most of the properties discussed below do not depend on the pre- 
cise structure of H, it can be considered as being H — R d . For the 
Laplace- Fourier transform, we need to consider its complexified space H 
(H = H + iH) equipped with the hermitian extension of the scalar prod- 
uct: 

VqE~H,\/xEH: (q,x) = (Re(q),x) + i(Im(q),x) (2) 

whereas for a Fourier transform it suffices to consider iH, and for a Laplace 
transform a subspace H + C H, e.g. H + = (R + ) d for H = R d . 

For any M G L(H,H), Syec( M) denotes its spectrum, Re(Spec(M) 
its real part, i.e. the set of the real parts of its eigenvalues, Ker(M) its 
kernel, sym(M) its symmetric part, i.e. 2.sym (M) = M + M* , where M * 
denotes the hermitian conjugate of M_. 



2 A generalized notion of scale 

The need for a generalized notion of scale will be more apparent in sectj^, 
however its definition is rather straightforward g|, at least for the linear 
case which will be sufficient for the following. However, the subsequent 
developments can be skipped in a first reading. 

Definition 1 An application \\ . \\ from a real Hilbert space H onto the 
set of positive real numbers R + is said to be a generalized notion of scale, 
associated to one-parameter (linear) dilation T , when it satisfies the fol- 
lowing: 

• nondegeneracy, i.e: 

||£||=0< = >x = (3) 

• linearity with the dilation parameter, i.e: 

Va; G #,VA G R + : || Z y x \\= A || x || (4) 

• Balls defined by this scale are strictly increasing, i.e.: 

VA, A' G R + ; A > A' : B\ 3 By (5) 
where the balls B\ defined by the dilation satisfy (due to eq.^|): 

B x =Z x {B^) = {x£H;\\x\\<\} (6) 
and have the following frontier: 

dB x = {x£H;\\ x ||= A} (7) 

It is straightforward to check that the canonical norm | . | is the 
scale associated to the isotropic dilation = AL Whereas the two first 
properties are rather identical to those of a norm, the last one is weaker 
than the triangular inequality which is required for a norm. 
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The conditions of existence of a generalized scale should depend on 
the generator G of the one parameter group of (linear) dilations : 

Z x = A= = e LnX = (8) 

Indeed, we have the following theorem B: 

Theorem 1 Let us consider the unit ball defined as an ellipsoid generated 
by a positive symmetric matrix A: 

Bi — {xe H; {x^.xf 12 < 1} (9) 

The dilation group T^ defines a generalized scale, if and only if its 
generator G satisfies: 

Spec{sym(A.G)) > (10) 
Indeed, the balls B\ are then defined by: 

A X = (T- 1 )*.A : .Z- 1 ■. B x = {xeH;(x,A x .x) 1/2 <l} (11) 
and we have: 

^fe4-£) = -(Sa 1 -^ sym{A.Q.TT x \x} (12) 

This theorem has the following corrolary, which will be sufficient for 
the following: 

Corollary 1 When the unit ball is an ellipsoid defined by a positive sym- 
metric matrix A ( eg.J^J, the dilation group defines a generalized scale, 
if and only if its generator G satisfies: 

Spec(sym(0) > Re{Spe(Q) > (13) 
and when A belongs to a given neighbourhood of a scalar linear application, 
i.e. 3fi € R + : A — fj.1. 

Indeed, for the scalar case (A = fiV), it is a straightforward conse- 
quence of theo.^. Due to continuity, it is also valid for a given neigh- 
bourhood of this case, i.e. when the eigenvalues of A are not too much 
different. However, simple counter examples (e.g. in L(R 2 , R 2 )) are easy 
to find, i.e. generators G which satisfy eqjlla, but violate eqltj as soon 
as the eigenvalues of A are different enough Jfj|. This shows that this 
neighbourhood is indeed bounded. 



3 Levy stable vectors 

Let us recall that a (real) Levy stable variable [], |^, [i| can be defined 
in the following manner, e.g. M (= d denotes the equality in distribution): 

Definition 2 A random variable X is said to be a Levy stable variable, 
iff it is stable under renormalized sum -i. e. with the rescaling factor a(n) 
and recentring term y(n) - of any n of its independent realizations Xi,(i — 
l,n). This corresponds to: 



3 



X l = d X(i = l,n) : Vn G N,3a(n),j(n) G R^X l = d a(n)X + 7 (n) 

! = 1 

(14) 

Furthermore, X is said strictly stable |j when the recentring term "i(n) is 
0. 

It is straightforward to check (by induction) that this definition is 
equivalent to the original definition given by Levy jl| which addresses the 
stability under any linear combination: 

Definition 3 Two identically distributed random variables X\,X 2 are 
said to be stable under linear combination iff: 



Xi = d X 2 = d X : Vai,a 2 G R + ,3a £R + ,y€R: a 1 X 1 +a 2 X 2 = d aX+ 1 

(15) 

Furthermore, they are said strictly stable, when the recentring term 7 
is 0. 

Due to their linearity, these definitions can be extended in a rather 
straightforward manner to random vectors |^, &, with the only neces- 
sary modification that the (rather trivial) recentring term (j(n) in eq.[l4|, 
7 in eq.[l5|) is now a vector Q However, this extension seems rather restric- 
tive, one of its consequences, recalled below, is that all the components 
have the same Levy index. We will therefore consider the following more 
intrinsic vector definition of a stable vector: 

Definition 4 A random vector X_ valued on an Hilbert H is said to be 
a Levy stable vector, iff it is stable under renormalized sum -i. e. with a 
rescaling linear application gin) and recentring vector 7 (n)- of any n of 
its independent realizations X_.,i = l,n. This corresponds to: 



Ki = d X(i = l,n) : Vn e N,3g(n) G L{H,H),j{n) G H : a(n).X+7_(r 

j=i 

(16) 

Furthermore, X_ is said strictly stable, when the recentring term j(n) is 0. 
The classical definition ^, ^] corresponds to the scalar case for a(n): 

g(n) = a{n)l (17) 
and therefore can be called 'quasi-scalar case' of Levy stable vectors. 



1 according to Feller's terminology 

2 In factjthis definition might have seemed so trivial to Levy that it is not explicitely 
written in |p| 
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4 Attractivity 

We can now extend the definition of attractivity to random vectors : 

Definition 5 The distribution R of the independent random vectors X f 
belongs to the domain of attraction of a distribution R a of the random 
vector X_, if the distribution of their renormalized sum -i.e. with a rescaling 
linear application a(n) and recentring vector b(n)- tends to R a . This 
corresponds to: 

7- L~d%= In — 1 —\/ d /io\ 

hm n ^ x — = X (18) 

ayn) 

By its very definition (def^J) each stable vector belongs to its own 
domain of attraction and we will show there is no other possible limit. This 
will demonstrate that the stable vectors are the only attractive vectors. 

5 Scaling law of the rescaling factor 

As for the scalar case, it is straightforward to demonstrate the following: 
Lemma 1 the rescaling factor a(n) forms a multiplicative group: 

Vm,n £ N : a(m.n) = a(m)a(n) (19) 

In order to obtain this group property, it suffices to iterate Eq.^| over 
m.n since it yields: 

a(rn.n)X_ + 7(m.n) = a(n).a(m)X + a(n).'y(m) + rwyln) (20) 
Lemma |l| yields the following: 
Lemma 2 The multiplicative group a(n) has a generator a -1 , i.e.: 

3g 6 L(H, H),Vn G N : g(n) = n=^ = e Ln[ - n ^ X ( 2 l) 
Due its symmetry in m and n, Eqj20] also shows that: 

36 eH,Vn: <y(n) = \a(n) - n] b (22) 
therefore, it demonstrates the following 

Lemma 3 If X f are stable with generator a~ and 1 ^ Spec(gt), there 
exists b ( eq$23[), such that 2Li~k are strictly stable with the same generator. 

6 Characteristic functions for stable vec- 
tors 

The first (resp. second) characteristic functions Zx_(q) (resp. Kx_{q)) are 
defined in the following way: 

Zx(q) = e K ^ = [ e^dPx (23) 
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where Px_ is the probability of X_ and the domain H' of the conjugate vec- 
tor q is respectively H' = H + ,iH,H for Laplace characteristic functions, 
for Fourier characteristic functions, and Fourier-Laplace functions^]. 

With the help of the lem.^, we need only to consider the case of strictly 
stable vectors. In this case the group property of a(n) (eq.^l|) corresponds 
to: 

Theorem 2 The second characteristic function Kx (<?) of a strictly stable 
stable vector has the following scaling behavior 

\/XeR + ,\/qeH' : Kx{Z x •?) = A Kx (g) (24) 

where the dilation T A has the generator cT 1 . 

For any positive integer A, this is an immediate consequence of def.^j 
and lem.[| It is readily extended to any inverse of integer A, by considering 
the intermediate vector q' — T_ x _ 1 -q, and therefore to any rational A. 
Finally, due to the continuity of Kx_ (q), this is true for any positive real 
A. In the classical case, we have: 

a = «I; \\y\\=\y\ a (25) 

and on the other hand the component (u, X) of a stable vector (X_) along 
any given direction u ([ u \= 1) is a stable variable. Therefore, the char- 
acteristic of the former is obtained from the one of the latter, with the 
help of a positive measure dE'(it) of the directions u. With the help of 
the lemj^, this yields the classical result Q 0, ^]: 

Corollary 2 The second characteristic function of a classical ( or quasi- 
scalar) Levy stable vector corresponds to: 

Kx(q)= [ (q,u) a dZ'(u) + (q,b) (26) 

where dS' is a positive measure which support dB[ is a subset of the 
frontier of the unit ball dB\ ^. 

it is straightforward that it is solution of egji^j. 

Eq.^ already points out a major difficulty with Levy stable vectors: 
a classical Levy multivariate distribution is in general not parametric, 
contrary to an attempt to reduce it to a parametric distribution |j] , which 
turns out [Q] to be only a very particular sub-case. 

One may note that in case of a Fourier characteristic function, eqj26| 
yields the classical expression for a £]0, 1[U]1, 2[: 



q = iq',q'£H; Kx(q) = i(j_,b) - {q,u) a dSf (u) + i0{j) (27) 

3 As for the scalar case the Fourier transform is defined for any type of stable vectors, 
whereas the Laplace transform is only defined for extremely asymmetric cases, but are more 
convenient for the latter case. 

4 The symmetry of the probability distribution is related to the one of dS, in the extreme 
asymmetric case, the support dB^ of this measure is a subset of H + fl dB\ . 
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where the asymmetry function (3(q) is given by: 



/3(g0 = tan(-) / (£u) | T" 1 dS'(u) (28) 

7 Levy canonical measure and generation 
of Levy stable vectors 

Let us recall that the Levy 'canonical' measure of a stable random variable 
can be best understood as corresponding to the distribution of hyperbolic 
jumps in a Poisson compound process Q [loj . 

This corresponds to substituting a Poisson sum for the deterministic 
sum in the definition of stability (eq.[l4| for stable variables, respectively 
eq jlq, for stable vectors), as well as for the definition of acttractivity 
(eq.|18|). The second characteristic function Kx_(q) is therefore rather eas- 
ily determined with the help of the is the Levy canonical measure which 
is the a— probability Fy_ of the jumps Y_. 

Indeed, a compound Poisson process of a random vector field Y_ defined 
by a measure da, generates the following (random) measures S_(A) for any 
given borelian A (of the space on which the process takes place): 



S{A) = / Yda (29) 

J A 

and which have the following characteristic function: 

Ks(A){q) = ca(A) [ [e^ y -> - u{(q,y)) - l]dF y (30) 



where the last term under the integral of the right hand side, is classical 
and merely removes the divergence of the a— probability F y in the limit 
| y [— » 0. Indeed, a <r — probability does need to be finite, but only to be 
the limit e — » of finite positive measures Fy. The term ui(q,y) rather 
corresponds to a generalization of the 'Levy's trick' for removing other 
divergences ((q,y) — > 0) of higher order (i.e. (q,y)), and whose appro- 
priate choice will be discussed below. In order to determine the latter, we 
have first to determine the scaling behavior of the Levy canonical measure 
Fy_. With no surprise, we obtain by conjugation of theorem |^: 

Theorem 3 The Levy canonical measure Fy of the jumps Y_ of a stable 
vector X_ has the following scaling behavior: 

VA G R + ,\/yeH : dFT'( y ) = \' X dFy_ (31) 
where the dilation T* is the conjugate °fH x , i- e - has for generator: a* -1 

In order to make a step further, let us consider the following generalized 
definition of the 'unitary vector' u(y) corresponding to a given vector y: 

Definition 6 In the framework of a generalized scale (deffy, the unitary 
vector u(y) of any non-zero vector y is defined by: 
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Vy e H - {0}, 3u(y) e 9Bi : = ^ (y) (32) 

In the case of the norm | . |, it corresponds to the usual notion of 
unitary vector, i.e.: 

u(y) = ||| (33) 

Theorem 4 The Levy canonical measure factors into a measure having 
a density in scale and a positive measure dT, on the frontier of unit ball 
dB x : 

d\\y\\ 

dF y _ = dnu)j^ (34) 

Indeed, due to its scaling property (theorem the Levy canonical 
measure Fy_ should factor into a scaling measure of || Y_ | (therefore a 
power law) and a measure invariant for any dilation T* . Since every ball 
B\ is obtained by dilation of the unit ball Si, it suffices to consider a 
measure on the frontier of the latter and to take care of the fact that 
the Levy canonical measure should be positive. The classical case (eq.p5[) 
does correspond to: 

dF y = Q^r y \~ a d Z(u) (35) 

\y\ 

and with the (classical) choice of either u>(x) = 0, (0 < a < 1) or uj(x) = 
x, (1 < a < 2), this yields eq j2q a nd the proportionality between the 
measure d£' (eqji|) and cffi (eq~||) is determined with the help of the 
Euler r function: 

dE , = rj3-a^ ds 
a(a — 1) 

As mentioned by Levy, another possibility is to choose ^| 

"M-ify (^) 

which is valid for any a G]0, 2], since ui ~ x(x — ► 0) and u) ~ 0(x — > oo). 

For the general case, the characteristic function is rather more in- 
volved: 

Corollary 3 The second characteristic function of a Levy stable vector 
has the following expression 



Kx{g) = / dE(H) r §[e ( ^l^ - Lo{q,r x {u)) - 1] (38) 
J ugbbx Jo 

5 One may note that Levy used in fact a simplified expression which is not relevant for the 
following 
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where lu is defined by egpfy dE is a positive measure on the frontier of the 
unit ball dB\ and the spectrum of a, which is the inverse of the generator 
of the dilation T , should satisfy: 



Eq.|3£j is a rather straightforward extension of the scalar case constraint 
(0 < a2). On the one hand, the lower spectral bound of a allows its 
inverse to be defined, whereas the upper spectral bound 2, which not 
surprisingly corresponds to the gaussian limit case, allows the definition 
of a generalized notion of scale associated to the generator a -1 (sect.^ and 
in particular theo.^) and it is required in order to ensure the convergence 
of the second integration involved in eq.^8| for the lower bound (A — ► 0). 
More precisely the integrand is of order: 



and it should be bounded by below by ^ to avoid divergences. We have 
the following two cases: 

• diaqonalizable a: 

eq.B9l merely corresponds to imposing the adequate constraint to 
each eigenspace, therefore to the full space. 

• otherwise: 

one needs to consider the Jordan form of a" 1 . And on each general- 
ized eigenspace corresponding to the eigenvalue a -1 = a+ib, a,b £ R 
the component (e.g.[jllj) of y — T_*^(u) is a linear combination of 

functions of the form | Log(X) \ k X a cos(bX) and | Log(X) \ k X a sin(bX), 
where k is bounded above by the codimension of the deficiency of 
the eigenvalue: k < c = inf ne M[codim(Ker((T_^ — Aj_) n )] Therefore 
the previous result holds. 

8 Conclusion 

We demonstrated that the notion of Levy stable vectors can be broadly 
generalized, the stability index becoming a linear application a which 
needs to respect only one spectral constraint (eq.^). We have two main 
cases: 

• diagonalizable q: 

on each of eigenspace, the stable vector have a common stability 
index which is the corresponding eigenvalue. However, there is no 
need for the eigenvalues to be equal, as in fact hypothesized in the 
classical definition of Levy stable vectors. 

• the generalization is in fact much broader: 

indeed complex eigenvalues of a induce rotation modulations, and 
deficiency of its eigenvalues introduce logarithmic modulations. As 
a consequence the stability index of a component of a stable vector 
depends on its norm. 



Re(Spec(a)) C]0,2] 



(39) 




(40) 
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Not only these results are fundamentaly important for multivariate 
analysis of random vectors, but for their stochastic simulations. Indeed, 
this study was in fact motivated by the latter, more precisely by simula- 
tions of multifractal fields jL| generated by strongly anisotropic stochastic 
differential equations. The latter naturaly introduce the generalized scale 
notion used in the present paper. 
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